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Electronic Properties and Persistent Spin Currents of Nanospring 

under Static Magnetic Field 

Taichi Kosugi 

Nanosystem Research Institute (NRI) "RICS", AIST, 1-1-1 Umezono, Tsukuba 305-8568, Japan 

Relativistic electronic properties of a nanospring under a static magnetic field are theoreti- 
cally investigated in the present study. The wave equation accounting for the spin-orbit interac- 
tion is derived for the nanospring as a special case of the Pauli equation for a spin- 1/2 particle 
confined to a curved surface under an electromagnetic field. We define the helical momentum 
operator and show that it commutes with the Hamiltonian owing to the helical geometry of 
the nanospring. The energy eigenstates are hence also the eigenstates of the helical momen- 
tum. We solve the equation numerically to obtain the surface wave functions and the energy 
spectra. The electronic properties are systematically examined by varying the parameters that 
characterize the system. It is demonstrated that either the nonzero spin-orbit interaction or the 
nonzero external magnetic field suffices for the occurrence of the persistent spin current on the 
nanospring. Two different mechanisms are shown to generate the persistent spin current. One 
employs the spin-orbit interaction coming from the local inversion asymmetry on the surface, 
while the other employs the curvature coupling with the external magnetic field. 

KEYWORDS: spin current, spin-orbit interaction, nanospring 
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1. Introduction 

Helices have been fascinating scientists for centuries 
not only in physics but also in other diverse fields since 
helical structures are widely seen in nature. One of the 
most famous examples is DNA, which has a double-helix 
structure formed by hydrogen bonds. Another example 
is proteins, many of which contain helical substructures. 
The helical structures of DNA and proteins are believed 
to play important fundamental roles in modern biology. 

On the other hand, the recent development of nan- 
otechnology facilitates the fabrication of materials in var- 
ious shapes and sizes, even including inorganic materials 
in helical forms. 1 They have various constituent com- 
pounds such as ZnO, 2-5 Si0 2 , 6-8 Pd, 9 SiC, 10 PbSe, 11 
and InGaN. 12 The length scales of such inorganic he- 
lices range from nanometers to micrometers. Since the 
helical form of nanomaterial allows it to be stretched 
or compressed without plastic elasticity if the force is 
not too strong, they are of interest for mechanical ap- 
plications. 13 The much more important and interesting 
properties for both theoretical and experimental stud- 
ies are the electronic properties, 14 ' 15 in which the helical 
and curved geometry should influence the behavior of 
electrons. Nanostructures with curved geometry in other 
than helical forms 16-21 have also been fabricated. 

Curved geometry is known to induce a geometric po- 
tential, 22 which affects the dynamics of an electron mov- 
ing on the curved surface, even when an electrostatic 
potential is absent. The curvature effects of surfaces in a 
helical geometry on the nonrelativistic electronic prop- 
erties without the spin degree of freedom have been 
theoretically studied. 23,24 There are many studies on 
examining the electronic properties of two-dimensional 
systems in other nontrivial geometries, which include a 
plane with a bump, 25 a cylinder in a transverse mag- 



netic field, 26 a catenoid, 27 and a rolled-up nanotube. 28 
Grigor'kin and Dunaevskii examined the electronic and 
optical properties of a cylinder with a helical potential in 
their works. 29-31 For one-dimensional systems, the elec- 
tronic and optical properties of a helix have been stud- 
ied. 32-36 The torsion-induced persistent charge currents 
on twisted quantum wires 37-40 and knotted tori 41 have 
also been studied. Despite the theoretical studies of the 
various kinds of nontrivial geometries reviewed above, 
no systematic study of the curvature effects of helical ge- 
ometry on the relativistic electronic properties with spin 
degree of freedom has been reported. 

The Schrodinger equations for a quantum mechani- 
cal particle on a curved surface have been used as tools 
for theoretical investigations on low-dimensional systems 
of nontrivial geometry. One of the most reliable meth- 
ods on which their formulations are based is the thin- 
layer method, proposed by da Costa. 22 This method re- 
gards the curved surface as a two-dimensional system 
embedded in a flat three-dimensional space. Ferrari and 
Cuoghi 42 adopted this approach and rigorously demon- 
strated, by choosing a proper gauge, that the separation 
of the on-surface and transverse dynamics under an elec- 
tromagnetic field is possible without approximation. Ko- 
sugi 43 extended recently their Schrodinger equation to 
the Pauli equation, which can describe a charged spin- 
1/2 particle with a nonzero mass confined to a curved 
surface under an electromagnetic field. 

In the present study, the relativistic electronic proper- 
ties of a nanospring 44 under a static magnetic field are 
examined. The wave equation accounting for the curva- 
ture effects and the spin-orbit interaction (SOI) is derived 
for the nanospring as a special case of the Pauli equation 
for curved surfaces. The equation is solved numerically 
and the interesting phenomena occurring on the spring 
are analyzed in detail. 
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This paper is organized as follows. In §2, the descrip- 
tion of the geometry of a nanospring as a curved surface 
is provided by using the curvilinear coordinates. In §3, 
the Pauli equation for the nanospring to be satisfied by 
the two-component wave functions is derived. In §4, the 
electronic properties of the nanospring arc systematically 
examined by varying the parameters that characterize 
the system. In §5, the conclusions of the present study 
are provided. 

2. Geometry of Nanospring 

2.1 Curvilinear coordinate system on nanospring 

In this subsection, the expression of the curved surface 
representing a nanospring is provided. Let us consider 
a circle of radius a on the xz plane whose center is at 
(R, 0, 0). By rotating the circle around the z axis while 
translating it in the z direction, we obtain a curved sur- 
face swept by the circle (see Fig. 1). Hereafter, we call 
it the nanospring, which is the two-dimensional system 
to be investigated quantum mechanically in the present 
study. An arbitrary point on the nanospring is repre- 
sented by the two coordinates 9 and cj> as 



(R + a sin (9) cos- 
(R + asind) sin< 
a cos 9 + T~f(f> 



(1) 



where — oo < <j) < oo is the azimuthal angle and 9 spec- 
ifies the position on the circle. We have defined 7 = 
where p > is the pitch of the nanospring. r specifics 
the chirality of the nanospring, which can take only +1 
or —I. 

The tangent vectors and the normal vector are given 

by 

/ cos 9 cos < 
or s ■ 

e e = ~K7T = a cos sin < 

80 \ - sin0 

dr s H*** 
— — — = 7 £ cos <p 



1 



e x 



£ sin 9 cos cf> + r cos 9 sin < 
£ sin # sin — t cos 6* cos | , (2) 
£ cos 9 



where £ 



_ -R+a sin /3 

7 



and rj = yj cos 2 9 + £ 2 . It is noted 



that eg and are not necessarily orthogonal. The sur- 
face metric tensor on the nanospring, gij — • ej(i, j = 
9, 6), and its determinant are calculated as 



a ) 9e<t> 



-Tajsin9, = 7 2 (1 + £ 2 ), 



g = det = a 2 j 2 r) 2 . 



(3) 



Hereafter, we use the dimensionless parameters R = R/j 
and a = a/7. 

2.2 Quantities needed for Pauli equation for nanospring 
A curved geometry is known to induce the geometric 
potential 22 in general, which affects the dynamics of a 
particle confined to the curved surface, whether the dy- 
namics is relativistic or not. This potential allows the 
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Fig. 1. (Color online) (a) Circle of radius a on the xz plane for 
generating a nanospring, viewed along the y axis. Its center is 
located at (R, 0, 0). By rotating the circle around the z axis while 
translating it in the z direction, the nanospring with pitch p is 
obtained, (b) Two examples of nanosprings with chirality r = +1 
on the left and t = -Ion the right. Solid curves on their surfaces 
represent constant-0 and -<j> contours. 



dynamics to be distinct from that on a flat surface. The 
geometric potential for an electron with effective mass 
m on the nanospring is calculated (see Appendix A for 
details) as 



V s (6) = - 



1 



2ma 2 



{£(1 + £ 2 ) +a{2 cos 2 9 + £ 2 ) sin(9} 5 



4?y 6 



a(£ 3 sintf — a cos 4 



(4) 



independent of (f>. Vs takes extreme values on the outer 
rim {9 — tt/2) and the inner rim (9 = 3ir/2) of the 
nanospring. 

To derive the Pauli equation for the nanospring, wc 
have to calculate the h matrices, 43 which are responsible 
for the spin-dependent part of the equation. They are 
calculated as 



1 



h <t>n = _ 



I 

arj 

1 



£ cos# 
-ie^C 

T 

ie^i 
sin 9 



ie~ l,p C 
— £ cos 9 



-le 

— T 



e-^cosi 
sin# 



(3) 



where £(#) = rcos^ — i^s\x\9. is antisymmetric with 
respect to the superscripts. For details of their derivation, 
see Appendix B. 
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3. Pauli Equation for Nanospring 

As an extension of the Schrodingcr equation for a 
curved surface under an electromagnetic field derived by 
Ferrari and Cuoghi, 42 the Pauli equation for the curved 
surface was provided recently by the author. 43 We derive, 
in this section, the Pauli equation for the nanospring as 
a specific case. 

3.1 Pauli equation in ordinary three-dimensional space 
Here we briefly review the Pauli equation in the ordi- 
nary three-dimensional space. It was originally derived 
by expanding the Dirac equation for a spin-1/2 par- 
ticle under an electromagnetic field using the Foldy- 
Wouthuysen method, 45 which is for the upper two com- 
ponents ^ of the four-component spinor. In the present 
study, we use the charge — e and the gyromagnetic factor 
g of an electron on the nanospring. The Pauli Hamilto- 
nian, which neglects the mass term mc 2 and the terms 
on the orders higher than m~ 2 , takes the form 



P 2m 



V + gnsS-B 



[U- S x E + S x E -U]. 



(6) 



4m 2 c 2 

II = p + -A is the kinetic momentum operator and 
the magnetic field B = V x A is the rotation of the 
vector potential. S = er/2 is the spin operator and <r 
is the Pauli matrix. V is the electrostatic potential, hb 
is the Bohr magneton. We identify m not with the bare 
mass m e , but with the effective mass of an electron. The 
quantization axis of spin is taken to be parallel to the z 
axis throughout the present study. 

3.2 Pauli equation and helical momentum of 
nanospring 

The author derived 43 the Pauli equation for a curved 
surface under an electromagnetic field by starting from 
the Pauli equation in the ordinary three-dimensional 
space of the form eq. (6). In the present study, we as- 
sume the electrostatic potential V to contain only the 
confinement part, which does not explicitly appear in 
the Pauli equation for the curved surface. The only po- 
tential felt by an electron moving on the nanospring is 
thus the geometric potential Vs , eq. (4) . The electric field 
is set to be perpendicular to the nanospring with a uni- 
form strength, and the magnetic field is set to be static 
and along the z direction: E = Ee n and B = Be( z y 
The electric field breaks the inversion symmetry on the 
surface and will be the origin of the SOI taking place 
in the present system. The magnetic field is realized by 
adopting the vector potential in the symmetric gauge, 
A = (B/2)(— yer x \+xer y \), which has the following com- 
ponents with respect to the curvilinear coordinates, 

Bi 2 e 



An 



0, A 



= - T ^l, (7) 



2 ' " 2r) 

on the surface of the nanospring. Since A n is nonzero, for 
establishing the surface Pauli equation, it is necessary to 
perform a gauge transformation 42,43 so that A n and its 
normal derivative vanish on the surface. Such a gauge 



transformation is always possible and Ag and A$ are un- 
changed via the transformation. The surface wave func- 
tion (f>) satisfies the time- independent Pauli equation 
for the curved surface as Hx = £Xj where H is the sur- 
face Pauli Hamiltonian and e is the energy eigenvalue. 
H is decomposed into the spin-independent and spin- 
dependent parts as H — ifsch + -Hsp— reij which are given 
by 43 



H — flsch + ^sp-rcl 



(8) 



Sch 



1 

2m 



i+e 

a 2 rj 2 



2 sinO 
ajr] 2 



a 1 if 7 



cos 9 I „ 9 a . , n \ _ 
- T j I l + e - -£sin0 )d 4 



— t f 4 1 + £ 2 + - ^ + 2 cos 2 0) sinf 

c larf 7 



+ 



2ie S7C 2 
c 2a?7 2 



rsmOde + -3$ \ - 



7 



e 2 ff 2 7 2 £ 4 
c 2 4ry 2 



(9) 



g/j, B B ia R 

-Hsp-rol — ^ O(z) 

2 ' a?] 



sin 9 
-e 1 ^ cos( 



t -ie'^i 
i&^i -t 



(10) 



where = + i eg ^ e - is the covariant derivative. 
a R = 4^z c i is * ne Rashba coefficient. 46 Owing to the 
nontrivial metric of the surface, the normalization con- 
dition of the surface wave function should be set as 

y/gdfM^x^X = const. (11) 

3.3 Helical momentum and its eigenstates 

In the present case, the electron traveling on the 
nanospring is inevitably forced to move circularly around 
the axis of the spring. Thus its movement can convey an 
orbital angular momentum. Here we define the helical 
momentum operator as 



L z + S z + r^p z 



-i(xd y - yd x ) + — ir'yd. 



-idd 



^, (12) 



where we have used eq. (1) for obtaining the last equal- 
ity. It is easily confirmed that K z commutes with the 
Pauli Hamiltonian ([H,K Z ] = 0), and we can thus ob- 
tain simultaneous eigenstates of the energy and the he- 
lical momentum. We therefore put the solution of the 
Pauli equation in the form 

y(6 6)- ^ (^^MO)) (13) 
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where v is real and can take continuous values. \ v is 
an eigenstate of the helical momentum with the eigen- 
value v. K z Xu = v\v It is clear from the definition of 
K v that the electron's circular movement in the positive 
(negative) direction of <j) provides positive (negative) con- 
tribution to v regardless of the chirality r. We adopt the 
normalization condition such that one electron resides in 
the nanospring per twist and the normalization condi- 
tion, eq. (11), becomes 

,•271 



d0(|^| a + |VHf) = l- 



(14) 



Whereas \ v has the dimension of inverse length, tjj v is 
dimensionless. From eqs. (9), (10), and (13), the Pauli 
equation becomes the following one-variable differential 
equation for the two-component wave function if) v with 
dimensionless parameters: 



Hv^v = e v ip v . 
The //-fixed Hamiltonian is given by 



(15) 



H u = - 



d9 2 



V 



(0) 

iy-l/2 





i+e d 2 
if 



Co) 



V 



(i) 



u-l/2 





V 



(1) 



V 



(1(9 



(16) 



where 



~ n\ sin 29 , o 



a£ sin I 



2a(v + B 1 £ 2 ) 



(17) 



Vv (0) = -^[2Ssin0£ 5 + {{a 2 + 4) cos 2 - 2}£ 4 

+25(1 -4 cos 2 0) sin6>£ 3 
+{(-6S 2 - 3) cos 4 9 + (55 2 + 9) cos 2 9 - 2}£ 2 
+4acos 2 0(cos 2 0-3) sin0£ - (2S 2 + 3) cos 4 9 + 5 cos 2 9] 
acos( 



+IT- 



2a£,sn\9{v - B 1 cos 2 9) 



-{v + B 1 ( 1 )(l + sin 2 9 + ( 2 ) 



9 



+ 2ma 2 V s , 



t —i£ 
i£ —t 



(18) 



a (-sm9{v-\ + B^ 2 ) cos9{v+\ + B^ 2 
i] [ cos 0{v - i + B^ 2 ) sin 9{v + \ + B 7 £ 2 ) 



W^ 1 



1 /r -i£ 



V \i£ — i 



We have defined the dimensionless parameters 

eB-f 2 - 



Q!r = 2maaB,, B~ 



2c 



-, £ z = ma 2 gfi B B. (21) 



It is noted that the eigenvalue e v is dimensionless and 
one can obtain the corresponding energy eigenvalue by 
multiplying e v by the energy unit e u = ^js- 

From eqs. (16) - (20), one can confirm that the Hamil- 
tonian for the helical momentum v under the magnetic 
field B and that with the opposite parameters —v and 
— B are related as a^ y )H^„{—B)*a^ = H U (B). Hence, 
if Vv(-B) is an energy eigenstate belonging to the eigen- 
value e v (B), criy^ip^B)* is an energy eigenstate for — v 
and —B belonging to the same eigenvalue: 

i>-u(-B) = <J iy) MBy,e-A-B) = e v {B). (22) 

Since a^ y )ip v (B)* is the time-reversed state of ip v {B), 
their spin directions arc opposite. The energy dispersion 
is even with respect to v when the magnetic field is ab- 
sent. 

One can also confirm that the Hamiltonians for 
nanosprings of opposite chiralities are related as 
H v {— t)* = H v (t), which implies 



Vv(-t) = 1p u (T)* , £ v (-t) = £ u (t). 



(23) 



3-4 Densities and currents of physical quantities on 
nanospring 

With the definition of the density and the current of 
an operator for a physical quantity, the plausible time 
development equation is obtained for the Pauli equation 
(see Appendix C). We define the density functions of 
the number, spin, and orbital angular momentum of an 
electron associated with the surface wave function Xv a s 

1 



2770777 



(hM' + hM 2 ). ( 24 ) 



Re(e^: t VM), 



= t O-(x) _ i_ 

S VX — Xv o Xv ~ r> 

I Zirajn 



s vy = xt G -f Xu = _i_l m (e^: t ^i) 



S VZ — Xv r> Xv — . 

2 47ra7?7 



(l^tl 2 " l^ll 2 ), (25) 



2ira 2 rf 



t sin 9lm('ij)^d9ip^ + VC4.de Vh) 




(26) 



respectively. While n v , s vz , and l vz do not depend on 
(p, s vx and s vy do. From the expression of s„, it is ob- 
vious that the spin direction at an arbitrary <f> on the 
nanospring can be obtained by rigidly rotating that at 
= with the same 9. It also means that s vx and s vy 
averaged over <f> vanish. 

Considering the definitions of currents, eqs. (C-7) and 
(CT3), we adopt the following definition of the current of 
an operator O associated with the surface wave function 
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30v — 3 Oh 



(27) 



jl v = —Re[ X lUO Xv + (H Xv yO Xv }, (28) 



; E x Re[xlaO X u] 



(29) 



4m 2 c 2 

Employing the relations between the energy eigen- 
states, eqs. (22) and (23), it is found that the densities 
and currents of the physical quantities for (r, B, is, </>), 
(r, — B, — is, (/)), and (— r, B, is, —<p) have the same magni- 
tudes. Their relative signatures are summarized in Table 
I. 

The </>-averaged total electron density is given by 

,-2k 



;£f) = ^ J J dvf{e ui ;e F )i 



(30) 



where e F is the Fermi level as a parameter and f{e;e F ) 
is the Fermi distribution function, i denotes the branch 
of the energy spectra. The </>-averaged total spin density 
s tot and the total orbital angular momentum density ^* ot 
are calculated similarly. s* ot and s^ ot obviously vanish. 
The total number of electrons as a function of the Fermi 
level is calculated as 

N tot (e F ) = [ d9^2nn tot {9;I F ). (31) 
Jo 

The total spin S tot and the total orbital angular momen- 
tum L* ot are calculated similarly. 

We define the net current of the j component of spin 
at a point as 



(32) 



The total spin current per pitch is calculated as its inte- 
gral over 4> and 9 as 

J t s f(e F )= [ y/gdO [ d<j>j s .(M;?f). (33) 
Jo Jo 



Table I. Relative signatures of densities and currents of physical 
quantities for (r, — B, —is, <f>) and (— r, B, is, —</>) with respect to 
those for (r, B, f, </>). 6 and the other parameters are assumed to 
be fixed. j cv represents the charge current. 



Quantity 



(t, -B, -v, 



{-T,B,V, 



III, 
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z + 

+ 

X — 
U + 
z — 
X — 
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z — 
x + 

V- 
z + 



4. Electronic Properties of Nanospring 

In the present study, the Pauli equation for the 
nanospring, eq. (15), is solved numerically. We adopt the 
periodic boundary condition with respect to 9, which is 
realized by expanding the two-component wave function 
in plane waves as 



With this expansion, 
(14), becomes 



E 



■inO 



(34) 



the normalization condition, eq. 



J2 i c «-ti 2 + i 



Cnv], I 



= l. 



(35) 



We take the terms only for |n| < 20 in the summation, 
which yields sufficiently converged results. The energy 
eigenvalues are thus obtained by diagonalizing the 82- 
dimensional (41 wave numbers for each of spin-up and 
-down states) complex matrix for each v. 

4-1 Energy spectra with no SOI and no magnetic field 

Let us first examine the energy spectra of the 
nanospring with the SOI and the magnetic field absent 
(a R = and B = 0). 

The dimensionless energies e„ of the nanospring of 
r = 1 for various combinations of p and R as functions 
of v are plotted in Fig. 2. The variations of the energy 
spectra due to the change in R/a look larger than those 
due to the change in p/a. Each of the energy eigenstates 
is found to be an eigenstate of <7( z ) unless the energy 
eigenvalue is degenerate. Even if the degeneracy occurs, 
linear combinations of the degenerate energy eigenstates 
for the eigenstates of <T( 2 ) are possible. This is reasonable 
since the only contribution that can mix the spin-up and 
-down components, the SOI, is absent in this case. 

In the case of large R/a, as shown in Fig. 2, the lowest 
two bands are close to each other and every four out of 
the higher bands appear in a group. This observation is 
understood by considering the limit of large R with a 
and p fixed, in which the spatial part of the eigenstate 
of H v is of the form e™ 9 with its eigenvalue proportional 
to n 2 apart from a constant, as seen in eq. (16). Hence 
the lowest group of two bands corresponds to n = and 
each of the other groups of four bands corresponds to 
±|n| 0. The spin degrees of freedom enter in addition 
and thus such groups of bands are formed. As seen in Fig. 
2, the energy dispersion for spin- up (-down) is symmetric 
around is = 1/2 (—1/2). The features of the dispersion 
will be discussed in detail in the next subsection by con- 
sidering the transformation laws of the Hamiltonian. 

The electron densities n v of the five lowest spin-up en- 
ergy eigenstates for various values of v are plotted in Fig. 
3(a), with the geometric parameters fixed at p/a = 2.5 
and R/a = 1.5. Among the eigenstates of a common 
branch, we observe the tendency that the electron den- 
sity near the axis of the nanospring for large is is small 
compared with that for smalle is. This is due to the de- 
pletion of the electrons from the inner region caused by 
centrifugal force. 
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p/a=2.5, R!a=\25 p/a=2.5, Rla=\.5 p/a=2.5, R/a=2.5 




-2-1012-2-1012-2-1012 




-2-1012-2-1012-2-1012 



V V V 



Fig. 2. (Color online) Dimensionless energy spectra of nanospring of r = 1 for various combinations of p and R as functions of helical 
momentum v. Magnetic field and SOI arc absent. Brighter red and darker blue curves represent purely spin-up and spin-down states, 
respectively. 
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The geometric potential Vs for various combinations 
of p and R is plotted in Fig. 3(b). It is seen that the 
influence of the change in R/a on the potential is much 
larger on the inner rim (8 = 3tt/2) of the nanospring 
than that on the outer rim. (8 = n/2) The larger R/a 
is, the weaker the 8 dependence of Vs is. This is easily 
understood by considering the asymptotic form of Vs- 
When R is much larger than p (R ^> 1), the geometric 
potential becomes 



V s (6) 



2ma 2 



1 

4 



2R 



(36) 



4-2 Electronic properties under static magnetic field 
with no SOI 

Hereafter, we set the geometric parameters of the 
nanospring as R = 45 nm, a = 30 nm, and p = 75 nm. In 
addition, we set the material parameters as m = 0.05m e 
and g = — 5, which are close to the values of bulk In- 
GaAs, 47 for simplicity of our analyses. 

Let us next examine the electronic properties of the 
nanospring under a static magnetic field with the SOI 
absent (ap, = and B / 0). The dimcnsionless ener- 
gies e„ of the nanospring of r = 1 under magnetic fields 
B = 1 and 2 T as functions of v are plotted in Fig. 4. 
Each of the energy eigenstates is an eigenstate of er^) 
as in the cases of B = above, since the SOI is not 
introduced. The band dispersion of spin-up (-down) is, 
however, no longer symmetric around v = l/2(— 1/2), 
which is due to the nonzero B. The introduction of the 
magnetic field not only caused the rigid shift of the band 
dispersion via the Zeeman term, but also the deforma- 
tion for the spin-up (-down) states through the terms in 



v„ _ 1/2 ana v u _ 1/2 



an d V^"jy 2 ) in the Hamilto- 
nian. It was found that the deformations of the dispersion 
from the B = cases to the B ^ cases in Fig. 4 are 
caused predominantly by the term a 2 (f + -B 7 £ 2 ) 2 /rj 2 in 

Vu°^ . This term acts as an effective potential consisting of 
three contributions: centrifugal {v 2 ), diamagnetic (B 2 ), 
and cross-term (vB) ones. Among them, the cross-term 
contribution is responsible for the asymmetric shape of 
the band dispersion of each spin. As seen in Fig. 4, the 
energy eigenvalues of spin-up (-down) states for positive 
v—1/2 (y + 1/2) are, on the whole, high compared with 
those for negative v — 1/2 (y + 1/2). This tendency is 
more obvious for B = 2 T than for B = 1 T. This can 
be understood by considering the classical dynamics of 
an electron moving at a velocity v on the nanospring 
under the magnetic field B, due to which the electron 
feels the Lorentz force —evxB. Thus a spin-up electron 
with positive (negative) v — 1/2, which corresponds to 
traveling in the positive (negative) direction of 4> regard- 
less of t, feels an attractive force toward (outward) the 
axis of the spring when B > 0, and acquires a higher 
(lower) energy due to the centrifugal potential. This is 
also the case for a spin-down electron with positive or 
negative v + 1/2. We should keep in mind, however, that 
this classical interpretation collapses when the spin cur- 
rent on the nanospring is considered, as will be discussed 
later. 



(0) 



It was confirmed that the asymmetric shapes of the 
band dispersion in Fig. 4 become their mirror images 
when the direction of B is reversed. On the other hand, 
their shapes do not change when only the chirality t is 
altered. These observations are consistent with the rela- 
tionships of eqs. (22) and (23). 



B=IT 




B=2T 



Fig. 4. (Color online) Dimcnsionless energy spectra of nanospring 
of t = 1 for R = 45 nm, a = 30 nm, and p = 75 nm under static 
magnetic fields. Material parameters are set as m = 0.05m e and 
g = —5. Brighter red and darker blue curves represent purely 
spin-up and spin-down states, respectively. 



Let us consider the electronic properties in detail from 
a mathematical viewpoint. The Hamiltonian, eq. (16), is 
now diagonal in spin space and the equation to be solved 
is decoupled for the spin-up and -down components: 



H v+ \/2 - Bz 




where 



H v = - 



\+e d 2 



T) 



de 2 



v de 



for which there is the relation 

H- u {-B)* = H V {B) 



(37) 



(38) 



(39) 



By using a one-component eigenfunction ^±1/2 an d its 
corresponding eigenvalue £^±1/2 of H v ±i/2, two simulta- 
neous two-component eigenstates of H v and CT( Z ) can be 
constructed as 



Vv : 



lpv-l/2 




V-Wl/2 



-1/2 



B 7 



e v = e 



v+1/2 



-B 7 



(40) 
(41) 
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Fig. 3. (Color online) (a) Electron densities n v of the five lowest spin-up energy eigenstates for various values of v. Magnetic field and 
SOI are absent. Geometric parameters are fixed at p/a = 2.5 and R/a = 1.5. (b) Geometric potential Vs for various combinations of 
p and R. 
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It is thus clear that the energy dispersion of spin-up 
states is obtained by shifting that of spin-down states 
rigidly by 1 in the v direction and 2Bz in the e v di- 
rection. Specifically, if ip v is a purely spin-up eigenstate 
for an arbitrary v with the eigenvalue e„, the spin- flipped 
state 0{x)i&v is an eigenstate for v— 1 with the eigenvalue 
£ v -\ = e v — 2Bz- Similarly, if ip v is a purely spin-down 
eigenstate for an arbitrary v with the eigenvalue e v , the 
spin-flipped state uu.\i\) v is an eigenstate for v + 1 with 

the eigenvalue e v +\ = e, „ + 2Bz- 




v=0.5 


Bf2T\ 






ffltlT\\ 




/ / S«0T\ \ 








Fig. 5. (Color online) Electron densities n„ of the lowest spin-up 
energy cigenstates for various values of magnetic field B. SOI 
is absent. Geometric parameters are set as r = 1, R = 45 nm, 
a = 30 nm, and p = 75 nm. Material parameters are set as 
m = 0.05m e and g = —5. 



the outer region of the spring and pour into the inner 
region as B increases. This is because, as the magnetic 
field becomes stronger, the diamagnetic potential, which 
is proportional to B 2 and pushes the electron toward the 
inner region, becomes stronger competing with the cen- 
trifugal potential, which pushes the electron toward the 
outer region. For the spin- up states of v = 1/2, the cen- 
trifugal potential is absent and thus the localization of 
the electron density in the inner region of the spring is 
more significant than those for v ^ 1/2. 

The total number of electrons 7V tot , the z component 
of total spin 5*°*, and the z component of total orbital 
angular momentum L* ot are plotted in Fig. 6(a) as func- 
tions of the Fermi level £f- The eigenstates for spin- up 
electrons are occupied for B > prior to those for spin- 
down states, since the spin-up states are energetically 
more favorable due to the Zeeman effect. The filling of 
the bands in such a manner is reflected also in the curves 
of Sl ot . Since the cyclotron movement of the electron in 
the negative <fi direction is energetically more favorable 
than in the positive one in this case owing to the Lorentz 
force, as discussed above, the negative L* ot is observed. 

The ^-averaged total number of electrons ?i tot , the z 
component of total spin s* ot , and the z component of 
total orbital angular momentum ^ ot are plotted in Fig. 
6(b) as functions of 9 for e F = 2.5 and 4 with B = 2 T. It 
is seen that n tot is larger around 9 = 3ir/2 than around 
9 = 7r/2, which means that the diamagnetic contribution 
dominates over the centrifugal contribution in this case. 
The contribution to L' ot from the outer region of the 
spring is much larger than that from the inner region 
since the larger distance between the axis and the outer 
region leads to a larger orbital angular momentum. 

Jj.,3 Electronic properties under static magnetic field 
with SOI 

Let us now examine the electronic properties of the 
nanospring under the static magnetic field with the SOI 
present. (a R ^ and B ^ 0) The dimensionless energies 
e v of the nanospring of r = 1 for various combinations of 
B and cur are plotted in Fig. 7. Some of the degeneracies 
in the nonrelativistic band structures are resolved by the 
SOI. For ojr 7^ 0, each of the energy eigenstates is not 
an eigenstate of o~r z \ since the SOI mixes the spin-up and 
spin-down components of the wave function [see eqs. (19) 
and (20)]. When the magnetic field is absent, the energy 
dispersion is even with respect to v also for or ^ [see 
eq. (22)]. 



The electron densities of the lowest spin-up energy 
eigenstates for various strengths of the magnetic field 
are plotted in Fig. 5. Although the spin- up electron den- 
sities for v = —1.5 and 2.5 are the same when B = 
since H*_ 2 = H 2 [see eq. (39)], when B / 0, they are 
different, as clearly seen in Fig. 5. The cyclotron move- 
ment of the electron with v = 2.5 is, as discussed above, 
enhanced by the magnetic field and thus the localization 
of the electron density in the inner region of the spring 
is stronger than that for v = —1.5. It is observed that, 
for all the plotted i>'s, the electron densities spread away 



4-4 Persistent spin currents 

From the viewpoint of applications for nanodevices in 
spintronics, we are interested particularly in the spin 
transport along the nanospring. We therefore examine 
here the behavior of the persistent spin current occur- 
ring on the nanospring in detail by varying the strengths 
of the SOI and the external magnetic field. 

The z components of the total spin current J* ot occur- 
ring on the nanospring are plotted in Fig. 8(a) for vari- 
ous combinations of or and B as functions of the Fermi 
level. These results indicate that either a nonzero SOI 
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Fig. 6. (Color online) (a) Total number of electrons N tot , z component of total spin S* ot , and z component of total orbital angular 
momentum L* ot as functions of Fermi level ep are plotted as solid curves for various magnetic fields B. Brighter red and darker blue 
dashed curves represent spin-up and spin-down contributions, respectively. The origins of energy are set to the lowest eigenvalues for 
the individual energy spectra. SOI is absent. Geometric parameters are set as r = 1, R = 45 nm, a = 30 nm, and p = 75 nm. Material 
parameters are set as m = 0.05m e and g = —5. (b) </>-averaged total number of electrons n tot , z component of total spin s* ot , and z 
component of total orbital angular momentum Z' ot are plotted as functions of 8 for ep = 2.5 and 4 with B = 2 T. 




Fig. 7. (Color online) Dimensionless energy spectra for various combinations of B and or as functions of v. Dashed curves are for 
apj, = 0. Geometric parameters are set as r = 1, R = 45 nm, a = 30 nm, and p = 75 nm. Material parameters are set as m = 0.05m e 
and g = —5. 
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or a nonzero magnetic field suffices for the occurrence 
of the persistent spin current. The z components of J* ot 
and J'*° t were found to vanish irrespective of the values 
of or and B. The x and y components of «7* ot also van- 
ish. (J t g)z oscillates in a complex manner around the 
origin as the filling is increased. This behavior clearly 
contradicts the interpretation introduced above for the 
asymmetric features of the band dispersion in Fig. 2. The 
classical picture employed in the interpretation led to the 
inequivalence of the positive and negative directions of 
4>. If this interpretation were also true for the spin cur- 
rent, the direction of ,/!j ot would not change even when 
the filling is varied. To capture the behavior of the spin 
current in detail, let us observe (J s ) z on the nanospring, 
as shown in Fig. 8(b). The regions of positive and nega- 
tive values exist irrespective of the sign of the total spin 
current. This result implies that the oscillatory behav- 
ior of the spin current on the nanospring is a quantum 
mechanical effect that is seen also in other systems. 48,49 

The mechanisms for the occurrence of the persistent 
spin current on the nanospring can be explained qualita- 
tively as follows. 

We first consider a situation in which the spin-orbit 
coupling is present and the external magnetic field is 
absent. In analogy with the Rashba effect 46 on a flat 
plane, two electrons traveling in the same direction with 
opposite spins are energetically inequivalent due to the 
SOI, as shown in Fig. 9(a). These two electrons thus give 
a nonzero contribution to the net spin current. This is 
also the case for two electrons traveling in the opposite 
direction, whose contribution is the same in magnitude 
and sign as that from the former two electrons. The four 
electrons depicted in Fig. 9(a) thus cause a nonzero spin 
current in total. 

We next consider a situation in which the spin-orbit 
coupling is absent and the external magnetic field is 
present. In this case, two electrons traveling in the same 
direction with opposite spins are energetically inequiv- 
alent due to the Zeeman effect, as shown in Fig. 9(b). 
These two electrons thus give a nonzero contribution to 
the net spin current. This is also the case for two elec- 
trons traveling in the opposite direction, and their con- 
tribution has the opposite sign to that from the former 
two electrons and has a different magnitude due to the 
inequivalence of the positive and negative <p directions. 
The four electrons depicted in Fig. 9(b) thus cause a 
nonzero spin current in total. 

The mechanism in the case of cur 7^ and B = 
explained above does not originate in the curvature of 
the nanospring. This mechanism is essentially the same 
as that for the persistent spin current on a flat sur- 
face. 50 The mechanism for qr = and B ^ 0, on the 
other hand, originates in the curvature of the nanospring, 
which couples with the external magnetic field to cause 
the inequivalence of the orbital motion in opposite direc- 
tions. 

(J t s ot )z for four combinations of the chirality and the 
direction of the external magnetic field are plotted in 
Fig. 10. It is seen that J* ot 's for the opposite r are in 
the opposite directions, while those for the opposite B 



(a) a R j±0,B=0 

z 




Fig. 9. (Color online) Schematic illustration of the mechanisms 
for the occurrence of persistent spin current on a nanospring. (a) 
When aji ^ and B = 0, two electrons traveling in the same 
direction with opposite spins are energetically inequivalent. They 
give a nonzero contribution to the net spin current. This is also 
the case for two electrons traveling in the opposite direction, 
whose contribution is the same in magnitude and sign as that 
from the former two electrons, (b) When a^i = and B ^ 0, 
two electrons traveling in the same direction with opposite spins 
are energetically inequivalent. They give a nonzero contribution 
to the net spin current. This is also the case for two electrons 
traveling in the opposite direction, whose contribution has the 
opposite sign to that from the former two electrons and has a 
different magnitude. 

are in the same direction. The magnitudes of J° were 
found to be the same for all the four combinations. This 
result is a direct consequence of the transformation laws 
of the wave function [see Table I and eqs. (22) and (23)] 
originating from the helical geometry of the nanospring. 
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Fig. 10. (Color online) z components of J* ot for four combina- 
tions of t and B as functions of Fermi level ep . Geometric param- 
eters are set as -R = 45 nm, a = 30 nm, and p = 75 nm. Material 
parameters are set as m = 0.05m e , g = —5, and or = 20 meV 
nm. 



5. Conclusions 

In the present study, we first derived the Pauli equa- 
tion for the nanospring to be satisfied by the two- 
component wave functions. The electronic properties of 



12 J. Phys. Soc. Jpn, 



Full Paper 



Author Name 



(a) 



(Ji°% (a.u.) 



_B=0T 


— «r = meV nm 

— «r = 5 meV nm 

— a R =10 meV nm 

— a R =20 meV nm 
















_B=17 











(b) ct R =20 meV nm, 5=0 T 
e F =0.25 



0.006 
0.004 
0.002 


-0.003 
-0.004 
-0.006 



-0.002 
-0.004 
■0.006 



1 



(c) «r=0 meV nm, B=l T 
e F =0.3 



0.003 
0.002 



-0.001 
-0.002 
■0.003 



e F =l.l 



%=0.75 



1 



s F =2.2 




Fig. 8. (Color online) (a) z components of J x ° x as functions of Fermi level ep for various combinations of an and B. The origins of 
energy are set to the lowest eigenvalues for the individual energy spectra. Geometric parameters are set as r = 1, R = 45 nm, a = 30 
nm, and p = 75 nm. Material parameters are set as m = 0.05m e and g = —5. (b) z components of j s on the surface of nanospring 
for < 6 < 2n. 
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the nanospring are then systematically examined by 
varying the parameters that characterize the system. 
The overall features of the band dispersion were demon- 
strated to admit the classical interpretation employing 
the Lorentz force acting on the electron. The spatial 
distribution of the electrons on the nanospring was in- 
terpreted to be a consequence of the competition be- 
tween the centrifugal and diamagnetic potentials. It was 
demonstrated that either a nonzero SOI or a nonzero 
external magnetic field suffices for the occurrence of 
the persistent spin current on the nanospring. Although 
we found that the behavior of the spin current on the 
nanospring does not to allow for the classical interpreta- 
tion, we were able to have the simple explanations of the 
two different mechanisms for the occurrence of the per- 
sistent spin current. One employs the SOI coming from 
the local inversion asymmetry on the surface, while the 
other employs the curvature coupling with the external 
magnetic field. A large part of the interesting phenomena 
observed in the present study comes from the helical ge- 
ometry, which specifically means that an electron moving 
circularly is inevitably forced to travel in the vertical di- 
rection. Similar effects on the electronic properties should 
thus be observed in other systems in helical or twisted 
geometries. The present work will evoke much interest in 
various curved systems for theoretical and experimental 
studies in the future. 
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Appendix A: Geometric potential of nanospring 

By taking derivatives of both sides of e„ • e„ = 1 with 
respect to q a (<ia = 0,<j>), we have -e n = 0. This means 
that the derivatives of e„ are linear combinations of the 
tangent vectors. Hence, from eq. (2), we can calculate the 
Weingarten matrix 22 a a b, which satisfies = a a bet, as 



1 



a cos 



cos 



> sin ( 



7 

£ sin 9 



1 



sm ( 



= t — , a<t,<h = 
rja 



(Al) 



By substituting a a b into the definition of the geometric 
potential, 22 



2m 



(Tra) 2 



det i 



(A2) 



that for the nanospring, eq. (4), is immediately obtained. 



Appendix B: h matrices of nanospring 

A point r close to the nanospring is represented by the 
three coordinates as 

r(6, <f>, n) = r s {6, 0) + ne n (6, 0), (B-l) 

where v is the distance between the point and the 
nanospring measured along the normal at (0,(j>). The 
derivatives of the Cartesian and the curvilinear coordi- 
nates with respect to each other, 

dxi„ , 

(B-2) 



% ' e w dx (a) ' 



satisfy the conditions = Sf, ef^e 1 ^ = 5^ due to 



the chain rule of derivative, ef^ on the nanospring can 
be calculated from eqs. (1) and (2), which then enables 



one to obtain 



» 



(«) 



as the inverse of the 3x3 matrix of 



as follows: 



'0*0 



(1 + £ 2 ) cos cos 4> — t£ sin 9 sin <j> 

2 ' 

(1 + £ 2 ) cos 9 sin + t£ sin 9 cos <fi 



e (z)="- 



£ 2 sin (9 



(x) " 
<P 

e (y) 



— £ sin 4> + t cos 9 sin 9 cos < 
IV 2 

£ cos 4> + t cos 9 sin 9 sin <fi 
j, cos 2 9 



IT 



£ sin 9 cos (f> + t cos 9 sin < 
V 

£ sin 9 sin <f> — r cos 9 cos < 
V 

£ cos 9 



e W = 



(B-3) 



By substituting these expressions into the definition of 
the h matrices, 43 h % i = e(abc) <T (o) e (; ) ) e '( c ) = where 
e (abc) is the Levi-Civita symbol, we obtain eq. (5). 

Appendix C: Definition of current of an opera- 
tor for generic Hamiltonian 

In this Appendix, we look for a possible definition of an 
arbitrary operator consistent with its time development 
equation. 

Here we consider a Hamiltonian for a two-component 
spinor ip of the form 



H-— H 

2m 



where H has the generic form 

H= E vAFf\r)-n + F?\r)]. 

j=0,x,y,z 

do is the identity matrix, ip obeys the Schrodinger equa- 



(CM) 
(C-2) 
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tion i^- = Hip. We assume F^ to be real. Let us con- 
sider the time development of the density function of a 
physical quantity represented by an operator 0, 



0(r,t) = Reft^Cty]. 



(C-3) 



We assume that O does not depend explicitly on t. The 
time derivative of this function is, from the Schrodinger 
equation, 

t 



.80 
dt 



tt2 / tt2 

tfO— i> - — </> ] O^j 
2m \ Im } 



We define the current associated with O as 



(C-4) 



(C-5) 



where we have defined the velocity operator as v 
—i[r,H]. From eqs. (CT) and (C-2), we have 



v = i+ E (c-6) 

Substitution of eq. (C-6) into eq. (C-5) leads to the de- 
composition of the current into two parts: j = +Joi 
where 

Jo = ^Re[V f nOV + {U^Oij] , (C-7) 



(C-8) 



The divergence of the kinetic part is calculated as 



p-jo = — imh^n 2 ©^ - (n 2 vo t £V] 

2m 



O. 



2m 



E (v 



F {1) +2lmF\ 0) )Re[^a J O^\, 



(C-ll) 



where O = — i\0, H]. The second and third terms on the 
right-hand side of this equation are the source terms. 
When they are absent, eq. (C-ll) is nothing but the con- 
tinuity equation. 

Let us consider here the expressions of currents for 
the Pauli Hamiltonian, eq. (6), as a special case of the 
definition, eq. (C-5). Comparing eqs. (6) and (C-2), we 
find 



(o) _ 



= 0,F 



6 e p 



V. 



Am^c 2 lJk 



-g^B r (C-12) 



We have used the Maxwell's equation Vx£ = — i O^f, 
which vanishes since the magnetic field is static in the 
present study. The third term on the right-hand side of 
eq. (C-ll) vanishes in this case. The residual part of the 
current hence takes the following form, from eq. (C-8), 



Jo = 



; E x Rc[^ trOip}. 



(C-13) 



4m 2 c 2 

For = 1, O = ip^ip is the electron density and its inte- 
gral over the entire space is the electron number. It is a 
conserved quantity since the source term vanishes and eq. 
(C-ll) reduces to the well-known continuity equation of 
the probability density. The integral of the density func- 
tion of an arbitrary operator is, however, not a conserved 
quantity in general. If O and H commute (O = 0), O is 
conserved. 

The definition of the spin current is mentioned here. 
Since the spin operator is the constant hermitian matrix, 
the spin current of the form of eq. (C-5) can be written 

as 



1 



(C-14) 



which is nothing but the conventional definition of the 
spin current. 



and that of the residual part j is calculated as 



pio = ilm 



+^a j Ff ) ■ pOi> + (p • FfWajOil) 



WW, 



zim 



J2(P ■ Ff ] ~ 2iImFj 0) )^a j Oip 



(C-10) 



By substituting eqs. (C-9) and (C-10) into eq. (C-4), we 
obtain 

80 
~dt 



-V-i + He[^6V] 
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